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Abstract 

Please note: after completion of this manuscript we learned that our 
main results, Theorem 12.11 and 12.21 can be obtained as a special case of 
Proposition 3.2 on Page 23 of Karlin's book [S]. 

Take a random variable X with distribution /x with some finite expo- 
nential moments, and weight it by an exponential factor e^'^ to get the 
distribution fi^ for the admissible 6'-values. Define also the so- weighted ex- 
pectation g{9) : — E^X with inverse function 6{g). This note proves that 
for a convex function -l-, E*('?'$(X) is a convex function of g, wherever it 
exists and is finite. Along the way we develop correlation inequalities for 
convex functions. Motivation for this result comes from equilibrium inves- 
tigations of some stochastic interacting systems with stationary product 
distributions. In particular, convexity of the hydrodynamic flux function 
follows in some cases. 
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1 Introduction 

Please note: after completion of this manuscript we learned that our main 
results, Theorem 12.11 and 12.21 can be obtained as a special case of Proposition 
3.2 on Page 23 of Karlin's book [B]. 
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Take a non-degenerate random variable X such that Ee < cxd for 9 £ I, for 
some open interval J. For these 9 define the exponentially weighted distribution 
'E%Y) : = E'^{Ye^^)/'Ee'^^ . The exponentially weighted expectation of X is 
g{9) : — E^X. This function is strictly increasing due to the nondegeneracy 
assumption. We denote its inverse by 9{g). 

Let $ be a convex function for which E^(^^$(X) exists in an open interval 
of £1- values. The first result of this note is the convexity of the function 

Motivation for this result comes from a class of asymmetric stochastic interacting 
systems that includes the zero range process (ZRP) and the bricklayer process 
(BLP). We explain these informally before turning to precise statements. 

The main application is related to the study of fluctuations of the current of 
particles J'-^-' (t) as seen by an observer moving at a fixed speed V. Equivalently, 
these are fluctuations of the height in the deposition formulation of the process. 
A key fact that underlies some of this work is that the variance of the current 
of a stationary process is linked to the deviations of a second class particle: 

(1.1) VariJ^^'Ht)) ^ C{g)E\Q{t) - [Vt]\. 

The variance on the left is taken in the stationary process at a fixed density g. 
The expectation on the right is taken with an altered initial product distribution: 
the density g invariant factor is put at each site other than the origin, while at 
the origin there is a different measure we denote by h'^^^\ defined in (|3.ip below. 
Q{t) is the position of a second class particle in the system. [Vt] is not the usual 
integer part but rather the integer between Vt and the origin that is closest to 
Vt. 

Identity (jl.ip has been known for the totally asymmetric simple exclusion 
process since the pioneering work of Ferrari and Fontes [5^. It was recently 
extended to the broader class of processes in [1]. For a complete discussion we 
refer the reader to [Ij, where the measures v^^^^l are defined in equation (2.6) 
and denoted by Jig. 

The presently relevant point is that for coupling purposes it is important that 
the measures v^^b) ^j-g stochastically monotone in the parameter g. In the case 
of the asymmetric simple exclusion process (ASEP) this is immediately obvious. 
This fact was utilized in a recent coupling-intensive proof ^ that established 
the order t^/"^ for the variance in (|1.1|) for ASEP when the observer travels at 
the characteristic speed. The first step in a program to extend the variance 
bounds to ZRP and BLP is to develop the coupling framework. However, the 
stochastic monotonicity of the measures jy^(e) is not at all obvious for ZRP and 
BLP. This we derive from the main convexity result in Section [3.21 

The conserved quantity in these asymmetric processes (typically viewed as 
particle counts, but also discrete gradients of the interface height) satisfies a 
hydrodynamic scaling limit where the limiting evolution is the entropy solution 
of a scalar conservation law of the form 

dtg + cd^H{g) = 0. 
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(We refer the reader to [7] for general theory.) The key quantity is the flux 
Ti.{g) which is computed as the expected jump rate in the stationary process at 
particle density g. The constant c is the mean increment of a particle that has 
decided to jump. 

A second application of the main convexity result is to give an alternative 
proof of the convexity of the hydrodynamic flux function 7i for zero range and 
bricklayer processes when the jump rate is convex. (See [1] for the original 
proof.) Concavity of the hydrodynamic flux also follows for concave jump rates 
in the zero range process. Strict convexity or concavity are also discussed. 

The characteristic speed referred to above is = cH'ig). So the issue can 
also be framed as the monotonicity of this quantity in the particle density. 

The rest of this note is organized as follows. We rewrite the convexity 
problem in terms of correlation inequalities of functions of X. These we handle 
via separation of positive and negative parts and further correlation inequalities. 
This is done in Section O In Section [3] we derive the consequences for stochastic 
interacting systems. 

2 Derivatives and correlations 

As in the introduction, let X be a nondegenerate real-valued random variable 
and $ a convex function on some interval that contains the range of X. The 
standing assumption throughout this section is that for some open interval / C 

R, 

(2.1) E(e^^)<cx) and E(X2|$(X)|e^^) < oo 

for all 9 £ I. Define the exponentially weighted distribution as E^(Y') = 
E^(re^'^)/Ee''-^. The function g{e) = E^'X is strictly increasing (justifica- 
tion below in CoroUarv 12. 4p . It has an inverse function 6{g) defined in some 
nontrivial open interval J. The expectation E^(^^<I>(X) is well defined for g £ J. 

Theorem 2.1. The function g t-^ E^(e)$(X) is convex on J. 

Our main interest lies in discrete distributions so we state a further condition 
for strict convexity for that case. Suppose the distribution ^ of X is supported 
on a discrete subset § of K. So S is either finite or countably infinite but locally 
finite. Let $ be a function defined on S. Extend $ to a function on the smallest 
closed interval that contains § by connecting adjacent points on the graph of 
$ with line segments. Assume the function $ thus defined is convex. Say $ 
is strictly convex at a point z G S if the slope of the extended <I> jumps at z. 
Such a point z cannot be the maximum or minimum of § because we have not 
defined 4> outside the smallest interval that contains S. If no such point exists 
then <& is linear. 

Theorem 2.2. Consider the discrete case described in the paragraph above. The 
function g t-^ E^'^^-'<i>(X) is strictly convex throughout its interval of definition 
if and only if S contains at least three points and $ is strictly convex at some 
point ofS. In the complementary case the function g ^ E^(^^$(X) is linear. 
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The remainder of this section covers the proofs. Throughout we only consider 
values 9 € I ioi which assumption (|2.1[) guarantees that the derivatives and other 
operations we perform are justified. In particular, since \X\'^ < fc!e~'^(e'^'^ + 
e~^'^), X has all moments under for each 6 E I because 6 zt e E I for small 
enough e > 0. 

We start with a preliminary lemma, repeated from Lemma A2 of [5]. 
Lemma 2.3. For any function tp, we have 

^EVW = Cov«(^(X),X) 

provided the expectations exist in a neighborhood of 9. 
Proof. 

— E^iX)-— ^ 

_Ei^{xyX^ E(X.e^^) 
^ B{^{X).e )• [Eg,x]2 

= Cov'{^iX),X). 



□ 

Recall that we exclude the degenerate case where /i is supported on a single 
point. 

Corollary 2.4. 

dg{9) _ d 



WX ^Cov^iX, X) ^Var^X >0, and 

UP uc 

(2.2) 



d9 d9 

dE«(e)^(X) dEV(^) d9{g) Cov''^<'\^{X), X) 



dg d9 dg Var«(^)x 

Now we proceed by rewriting the second derivative of E^'^)$(X) in terms of 
covariances. We omit the notation (X) from ^(X). 

Lemma 2.5. The following are equivalent: 

a) For any convex function E^^^^"I> is a convex function of g. 

b) For any convex function 

(2.3) Cov''($ -X, X)- Cov^(X, X) > Cov\<P, X) ■ Cov\X ■ X, X), 
where~ stands for centering w.r.t. E^. 



4 



Proof. We write, as in the corollary, 

dg Var^(«^'x 

We need to see if this is nondecreasing in p or, equivalently, nondecreasing in 9. 
That happens if and only if 

< Var^X ■ -^Cov^($, X) - Cov^($, X) ■ -^Var^X 
dp do 

= Var^X • [Cov^($X, X) - Cov''($, X) ■ E\X) - ■ Cov^X, X)] 

^ Cov''($, X) ■ [Cov^(X^ X) - 2E^X • Cov^(X, X)] 
= Cov^($X, X) ■ Cov^(X, X) - Cov^($, X) ■ Cov^iXX, X). □ 



Next we concentrate on proving that part b) of the last lemma holds for any 
distribution. Therefore we omit the superscript 9. 

Lemma 2.6. Part[6j] of Lemma \2.5\ is further equivalent to each of these two 
statements: 

c) For any convex function $ that is uncorrelated with X, Cov(<I>, X^) > 0. 

d) For any convex function $, 

(2.4) Cov($, X^) ■ Cov{X, X) > Cov($, X) ■ Cov{X'^, X). 

Proof. Given a convex function $, let ^{X) : — — C ■ X with C chosen so 

that $ is uncorrelated with X. $ is also convex, and we note that p.3p holds for 



4> if and only if it holds for $. Hence b) is equivalent to the statement obtained 
by restricting b) to convex functions that are uncorrelated with X. For such 
functions this statement becomes 

< Cov($X, X) - E$ • Cov(X, X) 
= Cov($, X^) + E$ • E{X^) - E($X) • EX - E$ • E{X^) + E$ • EX • EX 
= Cov($, X^). 

Thus b) is equivalent to c). 

Condition c) is a weakening of d), and we see that c) implies d) by deter- 
mining the constant in the transformation that led to <&: 

Cov($, X) = Cov($, X) - C • Cov(X, X) = 0, 

therefore 

^ Cov($, X) ^ 

$ = $ ^- L ■ X. 

Cov(X, X) 

Substituting this into Cov($, X^) > of c) leads to d). □ 
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Next we show that for part d) of the above lemma it suffices to consider the 
special case $(X) = 

Lemma 2.7. Part d) of the above lemma is implied by this statement: 
e) For any distribution (with finite third absolute moments) we have 
(2.5) Cov(|X|, X^) ■ Cov{X, X) > Cov{\X\, X) ■ Cov(X^ X). 

Proof. Consider functions of the form 

(2.6) (f>{x) = c + a ■ [x — xq]^ — b ■ [x — xq]^ 

for some a > b and xq, c G R. Notations + and ~ stand for positive and negative 
parts, respectively. These functions are convex. The first claim is that if (|2.4p 
holds for functions of this special form, then it holds for any convex 

This follows because $ can be approximated from below in a pointwise fash- 
ion by a sequence of functions of this type: 

m 

g{x) = c - ao[x - yi]^ + ai[x - + ^(a^ - afc_i)[a; - + 

fc=2 

with Og < ai < • • • < am and j/i < • • • < y™. The function g above is a 
sum of convex functions of type (|2.6p . To see the approximation, take points 
zq < zi < ■ ■ ■ < z„i and let be the slope of a tangent to <i> at the point 
{zi,^{zi)). Pick the z^'s so that the a^'s are strictly increasing. (This entails 
no loss of generality because a linear approximation to $ is exact throughout 
any interval with constant slope.) Let gi {0 < i < m) be the linear function of 
slope Oi that passes through the point {zi,^{zi)). Let yi {1 < i < m) be the 
x-coordinate of the point where the graphs of gi-i and gi intersect and set 

c = $(zo) + ao(yi - zq) = $(zi) + ai(yi - zi). 

Then it can be checked that g from above is the pointwise maximum of the g^'s, 
or equivalently, that g = gi on {yi,yi+i) with yg = -oo and y^+i = oo. By 
choosing the z^'s carefully one can create a sequence of convex functions g*-™^ 
such that g*^™^ y $ pointwise. By (|2.ip monotone convergence applies to show 
Cov(5('"), X'') -> Cov($, X'') for 6 = 1,2. 

Thus we can derive (|2.4[) for <i> by checking it for each g^™-*. Since (|2.4p is 
linear in $, it is then enough to know that it holds for each term of the type 
(j2.6p . This we now check. 

With suitably chosen constants A > 0, B and C, the transformation 

(2.7) (j){x) ^ A(l3{x) +Bx + C 

turns (f> of (|2.6[) into the function \x — xo\. (Note that a > 6 is needed for 
this.) The left and right-hand sides of (|2.4p are, up to the multiplying factor A, 
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invariant under these transformations. Hence (|2.4p holds for <j) if and only if it 
holds for |x — a;o|: 

Cov(|X - xo\, X^) ■ Cov(X, X) > Cov(|X - a;o|, X) ■ Cov{X^, X). 

Introduce now Y ^ X ^ xq, and write this inequality in the form 

Cov(|r|, (Y + xo)^) ■ Cov(y, Y) > Cov(|y|, Y) ■ Cov((r + xaf, Y). 

Subtracting 2xo ■ Cov(|y|, Y) ■ Cov(F, Y) from both sides leads to e) (for the 
distribution oi Y = X — xo) . □ 



Some elementary computations will now finish the proof of Theorem 12.11 
Lemma 2.8. Part e) in Lemma \27^ holds. 

Proof. For this proof, we introduce the positive and negative part moments: 
: = E((X+)% TV, :-E((X-)*). 

Expanding (|2.5p gives 

[P3 + - {Pi + Ni) ■ (P2 + N^)] ■ [P2 + N2- (Pi - Nif] 
> [P2 -N2- (Pi + A^i) ■ (Pi - A^i)] • [P3 -N3- (P2 + ^2) • (Pi - Ni)] . 

Somewhat tedious factoring shows that this is equivalent to 

(2.8) < TVi • (P3P1 - P2P2) 

(2.9) + Pi • {N3N1 - N2N2) 

(2.10) + P2N3 - P^P^N^ ~ P2N2N1 

(2.11) + P37V2 - P3N1N1 - P2P1N2. 

We proceed by showing that each line above is non-negative. Clearly if P{X > 
0} or P{X < 0} is zero, then P^'s or A^^'s are zero and the statement is triv- 
ially true. Assuming the contrary and dividing (|2.8p by [P{^ > 0}]^ makes 
conditional expectations out of the Pi's: 

^L-J^ = E{X^\X>0)-E{X\X>0)-E{X'\X>0)-E{X^\X>0). 

To show that this is non-negative, introduce the expectation 

p E{-xX^\X>0) 
E(X2|X>0) ' 

with which the previous formula becomes a constant multiple of 

EX ■ E— - 1 = -Covf X, — 

X \ X 
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Notice that the ^ measure is concentrated on positive values, where 1/X is a 
decreasing function of X hence the above covariance is non-positive. A similar 
argument shows that (|2.9p is non-negative. 
Separate (|2.10p into the sum of two terms: 

[P2P{X > 0} - Pi Pi] ■Nz + P2- [N^PiX < 0} - N2N1] . 

Divide the first bracket by \P{X > 0}]^ to get 

E(X2 I X > 0) - [E{X I X > 0)]2 > 0. 

Dividing the second bracket by [P{^ < 0}]^ leads to 

E(|X|3 \X <0)- E{X^ \X <0)- E{\X\ \ X <0) 



(2.12) 



= Cov(X^ \X\\X <0) >0 



since X^ is an increasing function of \X\ on non-positive numbers. The term 
(|2.1ip is treated in a similar manner. □ 

Tracing the lemmas backward shows that we have verified part a) of Lemma 
12.51 and thereby proved Theorem 12.11 

To prove Theorem l2.21 note first that in the complementary case ^{X) — aX 
on §, and then (|2.2p implies that the derivative dE^(^^$(X)/ dg is constant. 

To prove the main statement of Theorem 12.21 we retrace some earlier steps. 
Let G S be a point of strict convexity whose existence is assumed. Namely, 

(2.13) P{X <xo}F{X = xa}P{X>xo}>0 

and the slopes b — $'(xo — ) and a — ^'{xq+) satisfy a > b. Then we can write 

<^{x) = $(xo) + a{x - xo)+ - b{x - xq)^ + *(a;) 

for another convex function ^' that vanishes on an interval around xq. Since 
we already have Theorem 12.11 for 'J, it suffices to prove strict convexity of g t-^ 
E^(e)0(x) for 

(f>{x) = (f>(xo) + a{x - xo)^ — b{x - xo)~ ■ 



After an application of the transformation (j2.7p the question boils down to 
showing strict inequality in (j2.5p for the new variable Y = X — xq. For this 
it suffices to check that at least one of the quantities ()2.8p - (l2.1ip is strictly 
positive. From (j2.13p follows that each Pi and Ni is strictly positive. Schwarz 
inequality shows that P2P2 < P3P1 if Y has two distinct strictly positive values, 
and N2N2 < N3N1 if Y has two distinct strictly negative values. If both these 
requirements fail, then (|2.13p forces Y to take one positive value, one negative 
value, and the value zero with positive probability. But then this makes the 
quantity in (|2.12p strictly positive for Y. 

Thus we conclude that strict inequality holds in (|2.5p for Y = X — xq, and 
strict convexity of g 1-^ 'Ei^^'^^(f){X) follows. We have proved Theorem 12.21 



8 



3 Application to stochastic interacting systems 



To keep this note short we give a minimal possible introduction to the applica- 
tions of the convexity result and refer the reader to [1 and 4J for the complete 
picture. Let — oo < a;™'" < and 1 < x™^'^ < cx) be (possibly infinite valued) 
integers, and consider the discrete interval / = (a;™'" — 1, x™^^ + 1) n Z. Fix a 
function f : I ^ M+. For / 9 a:; > we set 

X 

/(^)!:=n/(y)' 

while for / 3 a; < let 

/(^)!:= ^ 



fi f{y) 

y=x+l 



finally /(O)! : = 1. Then we have 

/(x)!-/(x+l) = /(a; + l)! 

for all a; G /. Let 

log(liminf(/(^)!)i/^) 
oo , else 



a: = oo 



and 



log (limsup(/(-x)!)~^^'' ) , ifx™ 



-OO , else. 

We require / to be such that 9< < 6. In this case 



1 



5r f(yy- 

y£l 



defines a probability measure on /, and the exponentially weighted version 



is also well defined for any 9_< 6 < 9. This latter is the marginal of a stationary 
product distribution of many stochastic interacting systems, see e.g. [4]. 
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3.1 Convexity of hydrodynamic flux for zero range and 
bricklayer processes 

In particular, the attractive zero range process is an example where / = [0, oo)n 
Z, /(O) = < /(I), and / is non-decreasing. The rate for a particle to jump from 
a site with x particles is f{x). Its hydrodynamic (macroscopic) flux function 
H : K+ ^ M+ is given by 

W(^.)=E^(^)/(X) 

with the notation of the Introduction. The results of the previous section for / 
now read as follows: 

Proposition 3.1. // the jump rate f of the zero range process is convex (or 
concave), then the hydrodynamic flux H is also convex (or concave, respectively). 
Moreover, in this case H is strictly convex (or concave, respectively) if and only 
if f is not linear. 

The bricklayer process has / = (— oo, oo) n Z and / non-decreasing such 
that f{x) ■ /(I — x) = 1 for all x G Z. Its jump rate for a brick to be laid 
on a column between negative discrete gradients x on the left and y on the 
right is f{x) + f{-~y), see [¥ for more details. The hydrodynamic flux function 
H : K ^ M+ is now 

n{Q)^^'^'\f{X) + f{~Y)) 

where X and Y are i.i.d. variables with distribution Notice that non- 

decreasingness and non-negativity of f on Z excludes concave functions with 
the exception of the constant one function. Our result for this process is 

Proposition 3.2. // the function f of the bricklayer process is convex and not 
constant one, then its hydrodynamic flux H. is strictly convex. 

Parts of these two propositions were proved with coupling methods in [1] . 



3.2 Monotonicity of a special distribution 

We come to the primary motivation of the note. As explained in the Introduc- 
tion, the study of current fluctuations uses couplings of processes whose initial 
particle number at the origin obeys the following type of distribution: 

(3.1) v'^<^\y) = ] - E'^'^^] • /^'^(^) (2^™ < y < ^""n- 

Var ^^>X ^^^^ 

(See [H eqn. (2.6)] for the original definition.) To create couplings with useful 
monotonicity properties, one needs these distributions to be monotone in the 
parameter g, in the sense of stochastic domination. This we can now derive as 
a consequence of the main result. 
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Proposition 3.3. The family of measures v^^^) jg monotone in g. 
Proof. By Corollary [ 



^'(y) = ■ E([X - I^'iX)] ■ 1{X > y}) 

Var A 

^ Cov^(XJ^X>^ d 

Let us denote the i/^*^^^ -expectation by E'^'^^^^. Monotonicity of the family 

is equivalent to the property that, for any bounded non-decreasing function ip, 

< ^E''''^('?V(^)- 
d^i 

We compute a different expression for this derivative. Passing the derivative 
through the sum in the third equality below is justified because the series in- 
volved are dominated by certain geometric series, uniformly over 9 in small open 
neighborhoods. This follows from the definitions of ^ and and the assumption 
0<O<9. 



= E 'P{y)-^^[P'^'Hx>y}-l{0>y}] 

= -^E^(^) E ^(y) ■ > 2/ > 0} - 1{0 > y > X}] 



^E^(^) [E ^(y) - E ^(y)] - -^^'^'^my 

^ y=l y=X ^ 



X-1 



Above we introduced the function 

x-l 

$(a;) = E'^(y)-E^(y)' 

y=l y=x 

with the convention that empty sums are zero. To conclude the proof, notice 
that <i>(x + 1) — ^(x) = ^p{x). Thus a non-decreasing function ip determines a 
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(non-strictly) convex function <& with <&(!) = 0, and vice-versa. Hence Section 
[2] establishes that 

dg dg'' 
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